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Abstract 

The jump processes W(t) on [0, oo[ with transitions w — > aw at rate bw^ 
(0 < a < 1, b > 0, /3 > 0) are considered. Their moments are shown to 
decay not faster than algebraically for t — > oo, and an equilibrium probability 
density is found for a rescaled process U = (t + k)~@W. A corresponding 
birth process is discussed. 

PACS numbers: 
02.50.Ey Stochastic processes 
05.20.Dd Kinetic theory 
83.70.Fn Granular solids 
05.40.Fb Random walks and Levy flights 

AMS 1991 MSC : 
60C18 Self-similar stochastic processes 

82C23 Exactly solvable dynamic models (time-dependent statistical physics) 
70F35 Collisions 



1 



I. INTRODUCTION 



Random walks with absorbing states occur quite commonly. Here we focus on the asymp- 
totic behaviour of random walks suffering critical slowing down on approaching such a state 

HI- 

Consider a markovian random walk W on [0, +00 [ with continuous time t. The value of 
W is multiplied at random times by a given factor < a < 1, and successive jumps occur 
independently, with a waiting time between them distributed according to an exponential 
law with parameter bW 13 , where b > 0, j3 > 0. By introducing a new time t' = bt, a new 
variable W = W" and a new constant a' = a 13 , one reduces this model to the case (3=1, 
6=1 with no loss of generality. A characteristic parameter related to a is the e-folding 
number of jumps 

u = -l/\n(a) (1) 

such that, after cv jumps, W has been reduced by a factor a cv = e~ c . 

Denoting by f(w,t) the probability distribution function of W, the evolution equation 
for / reads 

d t f(w, t) = -wf(w, t) + a~ 2 wf{a~ l w, t) for a > (2) 
d t f(w,t) = -wf(w,t) + 6(w)x 1 (t) fora = (3) 
If they exist, the moments Xk(t) = (W k ) = J °° w k f(w,t)dw satisfy the hierarchy 

x k = -(1 -a k )x k+l (4) 

with the dot denoting derivation with respect to time t. We assume that P(W / (0) > 0) = 1. 

As this random walk can move only towards the origin, we are interested in its asymptotic 
approach to 0. Equations (Q)-(|4]) indicate that W and t _1 are dimensionally homogeneous, 
and no absolute time scale is available in these equations. 

Physically speaking, this random walk appears naturally in the kinetic theory of inelastic 
systems. Up to dimensional constants, if W is the modulus of the velocity of a particle suf- 
fering inelastic collisions (with restitution factor a) with fixed colliders, our process describes 
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the particle's velocity slowing down. For an assembly of such particles, the resulting process 
is thus inelastic cooling. However, further physical analysis of actual granular systems leads 
to a more complex process that our simple random walk [@,[3[]. 

II. EVOLUTION FOR W 

The linear hierarchy (|]) has a straightforward formal solution |l]U 

Xk{t) = 2^ -7j fT— jj, X*+l(0) (5) 

z=o ' 

with the q-factorial notation [A;]! = Ylt=i f(l ~ ° n ) /0- ~~ a )) ■ Note that lim a _,i[fe]! = k\ and 
that < limfc-Kx^l - a) k [k}\ < 1 for < a < 1. 

If Xfc(O) < M fc for some M > uniformly with respect to > 0, then (|^) implies 
that all moments Xk(t) are entire functions of t, bounded by M fc e ( - 1_a - )M ' t '. Also, if < 
f(w, 0) < cexp(— c'w 1 ) Ww g]0, oof for some c > 0, d > 0, 7 > 1, moments (fc e N) are 
entire functions of time. Conversely, hierarchy (^) is ill-defined if some initial moments are 
infinite, though (|2]) and are well-defined in this case too. These bounds are useful to 
discuss analytic properties of (H), but they are poor for positive times, as < < x^(0) 
obviously for t > 0. 

Consider now the evolution equation (Q). By linearity and scale invariance, its general 
solution reads 

f(w,t)= / f( Vj 0)Q(-,vt)- (6) 
Jo V V 

where Q is the solution of (0) for Dirac initial data Q(w, 0) = 8{w — 1). 

For a = 0, Q(w,t) = e _t <J(w-l) + (l-e-*)<J(tu), so that for 10 > 0, /0,t) = e^/^, 0). 
The moments follow readily, Xk{t) = k\t~ k ~ 1 f(0, 0) + o(t~ fc_1 ), and the surviving fraction of 
the initial population F(W(t) > 0) = i _1 /(0, 0) + o(t" 1 ), if /(«;, 0) is smooth near w = 0. 

For a > 0, the representation Q(w,t) = J^^Lo^^)^ 10 ~ a ™) yields coefficients q n by 
induction. Indeed, % = — Qo, Qn = —Oi n q n + a n_1 <?n-i, with initial data qo(0) = 1, ? n (0) = 
for n > 0. The Laplace transforms q n (s) = J °° q n (t)e~ st dt follow for 9R(s) > : 
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qo(s) = (1 + s)- 1 (7) 

- / x C^™" 1 - 1 TT Q;fc ~ 1 /ox 

s + a n s + l- LJ -s + a fe 

fc=i 

The accumulation of poles s n = —a n to the origin reflects the critical slowing down of the 
process. 

The Laplace transforms X k (s) of moments X k (t) = J2 n Qn(t)a nk have singular expansions 
for s — > 0. In particular, (1 + s)Xi(s) = 1 + Xi(s/a) for s > 0, which admits the solution 
Xi(s) = A(s)\ns + B(s) with entire functions A(s) = X]^Lo a ™ s ™' B{s) = J2^=o near 
s = : a n = a [ILi^ - * ~~ & « = ( & «-i ft " + a n lna)/(l - a n ), a = v. 

III. RESCALED PROCESS 

As the natural time scale for the evolution of W is 1/W, consider the rescaled variables 

r = ln(l + t//c) (9) 
U(r) = (t + K)W(t) (10) 

where the characteristic time k is adapted to the specific initial data /. The moments y k of 
U and its probability density h are 

y k (t) = (t + K) k x k (t) (11) 
/i(u,r) = (* + «)/(«;,*) (12) 
and for < a < 1 obey the equations 

Vk = ky k - (1 - a k )y k+1 (13) 

y fe (0) = z fe (0)« fe (14) 

d T h = Lh(u,r) + a~ 2 uh(a~ 1 u,r) (15) 

%,0) = k -1 /(u/«> ) (16) 

where Lh(u, r) = —ud u h(u, r) — (tt + l)/i(u, r) and the prime denotes derivation with respect 
to r. It is easily seen that U is a stationary Markov process, ergodic on the half line ]0, oo[. 



FIG. 1. For a = 0.8 (solid line), 0.7 (dots) and 0.6 (dash-dots) : (a) invariant density h eq vs 
u ; (b) moments vs k (in semi- log scale). 

It jumps down by a factor a with rate u and drifts upwards along exponential characteristics 
(it' = u). 

As [/ is ergodic, one finds a stationary density h ccl {u) and its moments y eq : 

«,^n o^^a -<*)-' (l7) 

Moreover, as (|T3D-(|I4|) hold for all k G R and yo = 1 by normalisation, one finds yj q = ^ = 
— 1/ ln(a) in the limit k —>■ 0. Induction for negative /c shows that the stationary solution of 
(|IB|) has also finite moments for all — oo < k < +oo (diverging for \k\ — > oo). 
A series for /z eq is found in the form 

oo 

h e \u) = ^r lm u' 1 e- u/am (18) 

m=0 

with coefficients r] m = —ri m -\or x {cr m — 1) _1 (note that r] m ~ a" 1 '" 1-1 " 2 for m — > oo). This 
series converges in the half plane 9?(w) > 0, with an essential singularity at u = 0. At that 
point, all derivatives of /i eq vanish, in agreement with the fmiteness of its moments. Figure 
[I] displays this stationary density and its moments. 

While initial distributions h(u, 0) relax to the stationary density h eq , the distribution 
of the original variable W approaches J/] this profile, up to the rescaling by (t + k). In 
particular, Figure Q displays the first two moments of the Green function Q and the leading 
term (t + K)~ k y c ^ in its asymptotic expansion for t — > oo. Two choices for time scale k 
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FIG. 2. Expectations X^{t) for the kernel Q with a = 0.8 and (a) k = 1, (b) k = 2. Solid 
line : direct sum of series (||) with Xfc(0) = 1 V/c ; thin lines : asymptotically leading approximation 
(t + K)~ k y cq with k = 1 (dots) and k = u (dash-dots). 

are compared : n = 1 may seem natural in view of initial data P(W / (0) = 1) = 1, but the 
choice k = yl q /Xi(0) = v has the advantage that the first moment satisfies yi(0) = y eq and, 
indeed, the relaxation of (t + n)Q(u/(t + n),t) to h eq is almost unnoticed on these lower 
moments with the choice k = v. 

IV. POPULATION DYNAMICS INTERPRETATION 

One may also interpret our process as a birth (only) process for the variable Z = 
ln(ftVl^)/ma (with the same characteristic time n as introduced in ([]) and (|TUD) : Z in- 
creases by unit jumps, with a jump rate K _1 a z '. This description suggests that Z would 
describe e.g. a self-inhibiting single-species population growth process, where all present 
individuals of the species are cooperating to produce one more individual. 

It is easily seen that, for the population corresponding to the asymptotic distribution, 

(Z(t)) w u\n(l +t/n) - z/(ln[/) eq (19) 

where (ln[/) cq = lim^o dy^/dk. The population grows to infinity, logarithmically in time, 
as does the solution z = uln(a + t/n) — v\nv (with integration constant a) to the rate 
equation z = K~ x a z corresponding to the birth process. 
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V. CONCLUDING REMARKS 



The long-time behaviour of moments in the case f(w,0) does not vanish on a neigh- 
bourhood of w = raises interesting questions, as the density h eq is flat at the origin. 
Relaxation of (t + K,)f(u/(t + n),t) is likely to be algebraic, and the competition between the 
scale factor (t + n)~ l and the relaxation of h may lead to non- universality in the asymptotics. 

One may also consider that the absorption at the origin is a coarse description of the 
relevant physical processes. In this respect, one may balance the drift towards zero by 
various simple processes : 

1. adding a diffusion term to (fj), with reflecting boundary condition (d w f(0,t) = 0) |§ ; 

2. imposing a constant acceleration (with W interpreted as a velocity) between jumps, 
as in the one- dimensional inelastic Lorentz gas model ||. 

Analytic investigations of these models reveal rich behaviours and varied structures of equi- 
librium distributions ||||. 
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